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industry,  where  it  is  essential  to  the  design  of  any  modern  commercial 
integrated  circuit,  research  on  compute! — aided  control  system  design  has 
been  held  back  by  the  primitive  electro-mechanical  and  hydraulic  hardware 
with  which  control  systems  have  historically  been  implemented.  Indeed, 
computer  aids  are  not  required  for  the  design  of  the  type  of  low  order 
(usually)  single  loop  control  systems  which  can  be  implemented  on  such 
hardware.  With  the  development  of  the  modern  microcomputer,  however,  it  is 
now  feasible  to  implement  a  sophisticated  high  order  multi-loop  control 
strategy  and,  as  such,  the  control  community  has  begun  to  look  towards  CAD 

with  renewed  interest.^ ’ 10 ’ 1 1 

Unlike  manual  design  techniques,  wherein  the  designer's  intuition  can 
compensate  for  the  lack  of  a  complete  analytic  theory,  a  CAD  package  must 
be  based  on  a  precise  analytic  theory  if  it  is  to  be  successfully 
implemented.  Indeed,  the  alternatives  of  brute  force  optimization  and/or 
disorganized  interactive  searches  are  neither  cost  effective  nor  "socially" 
acceptable  to  the  engineering  community.  As  such,  the  goal  of  the  research 
described  in  the  present  proposal  is  the  development  of  an  analytic  theory 
which  can  serve  as  the  basis  of  a  CAD  package.  In  particular,  we  are  in 
the  process  of  developing  a  complete  analytic  theory  for  the  incorporation 

of  the  fundamental  asymptotic  and,  hopefully,  initial  .value  constraints 

\ 

into  a  CAD  package.  Specific  research  topics  will  include: 
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(i)  pole  placement 

(ii)  initial  value  constraints 


(iii)  design  of  tine-varying  systems 

(iv)  simultaneous  design 

Our  theory  is  predicted  on  the,  now  classical,  stabilization  theory 
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developed  during  the  mid- seventies  by  Youla ,  et  al .  ’  ’  and  its 

extension  into  a  general  purpose  asymptotic  design  theory  as  formulated  by 

the  author,  Desoer,  Liu,  Vidyasagar,  Francis,  et  aj , 14 » ^ ,  1 1  , 1 2 , 1  3 , 1  ^  » 1 7 
Since  the  theory  is  now  well  known  we  will  simply  review  the  essentials 
here,  formulating  the  notation  which  will  be  required  in  the  sequel.  Foe  .. 
brevity  this  review  will  be  restricted  to  the  single  variate  case  though, 
in  fact,  most  of  the  theory  can  be  extended  to  the  general  linear  case. 

As  usual,  we  let  p(s)  denote  a  single  variate  rational  plant.  Here, 
p(s)  may  be  unstable  or  even  improper.  Of  course,  as  with  any  rational 
function  p(s)  can  be  represented  as  the  ratio  of  polynomials.  For  our 
purposes,  however,  we  find  it  more  convenient  to  express  p(s)  as  the  ratio 
of  stable  rational  functions,  n(s)  and  d(s). 


p(s)  *  n(s)/d(s) 


(1.1) 


As  with  the  classical  polynomial  representation  we  require  that  n(s)  and 
d(s^  be  copr i me ,  though  in  the  context  of  our  rational  fractional 
representation  it  suffices  to  only  require  that  n(s)  and  d(s)  h.;ve  no 
common  (closed)  right  half  plane  zeros.  Indeed,  "real  world"  control 
system  designers  often  permit  left  half-plane  pole /zero  cancellations. 


Interestingly,  this  coprimeness  condition  manifests  itself  in  the  Bezout 
equality 

u(s)n(s)  ♦  v(s)d(s)  =  1  (1.2) 

i.e.,  stable  rational  functions  n(s)  and  d(s)  are  coprime  in  the  sense  that 
they  have  no  common  right  half-plane  zeros  if  and  only  if  there  exists 
stable  rational  functions,  u(s)  and  v(s),  such  that  1.2  is  satisfied. 
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Given  the  above  notation  it  can  be  shown  ’  that  a  compensator,  c(s), 
will  stabilize  a  unity  gain  feedback  systems  with  plant,  p(s),.  as 
illustrated  in  Figure  1  if  and  only  if  it  takes  the  form 


-d(s)w(s)  +  u(s) 
.  "  'n(i)w(i)"‘v(i) 


(1 .3) 


where  w(s)  is  an  arbitrary  stable  rational  function  (which  serves  as  our 
design  parameter) . 

'■  -r 


the  form 


his)  «  -n(s)w(s)d(s)  +  n(s)u(s) 


(1.H) 


As  such,  one  can  formulate  a  design  theory  in  which  the  stability  question  # 

Is  totally  resolved  and  at  the  same  time  the  remaining  design  problem  is 
simplified  by  the  linearity  of  the  feedback  system  gain  expressions. 

Indeed  by  restricting  w(s)  to  stable  rational  functions  of  the  form  • 


w(s)  »  a(s)w(s)  +  b(s) 

where  a(s)  and  b(s)  are  appropriately  defined  stable  rational  functions  and 

w(s)  is  a  new  stable  rational  design  parameter  we  can  incorporate 
essentially  any  asymptotic  design  constraint  into  our  theory.  Specific 
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constraints  which  we  have  investigated  *  include: 

(i)  tracking 

(ii)  disturbance  rejection 

(iii)  robust  design 

(iv)  model  matching 

(v)  pole  placement 

(vi)  design  with  restricted  compensators 

(vii)  simultaneous  design 

The  goal  of  the  research  is  to  refine  the  above  described  results, 
especially  with  regard  to  the  pole  placement.  To  extend  the  theory  to 
permit  the  inclusion  of  initial  value  constraints,  and  to  investigate  the 
possibility  of  applying  the  theory  to  linear  time-varying  system,  via  the 
medium  of  a  recently  developed  generalized  frequency  domain  theory. 


i 


/> 


i 
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During  the  past  year  our  work  has  been  concentrated  on  the  solution  of 
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the  single  and  multivariate  pole  placement  problems^*  '  and  on  the 
extension  of  our  theory  to  the  case  of  linear  time-varying  systems. 

SINGLE-VARIATE  POLE  PLACEMENT:  By  the  pole  placement  problem^’1  **  in 
our  frequency  domain  setting  we  refer  to  the  problem  of  finding  a 
stabilizing  compensator,  c(s),  for  a  given  plant,  p(s),  such  that  the 
input/output  gain,  h(s),  for  the  feedback  system  of  Figure  1  takes  the  form 

h(s)  =  r(s)/q(s)  (2.1) 

where  q(s)  is  a  prescribed  Hurwitz  polynomial  and  r(s)  is  an  arbitrary 
polynomial  such  that 

o(r)  <  o(q)  (2.2) 

Here  o(  )  denotes  the  order  of  its  polynomial  argument  and  equation  2.2  is 
required  to  guarantee  that  h(s)  has  no  poles  at  infinity.  It  is 
interesting  to  compare  this  "frequency  domain  pole  placement"  problem  with 
the  classical  state-space  pole  (actually  eigenvalue)  placement  problem.  In 
the  latter  case  memoryless  state  feedback  is  employed,  thereby  fixing  the 
order  of  the  resultant  feedback  system,  whereas,  in  our  case  dynamic 
feedback  is  employed  allowing  the  designer  to  change  the  degree  of  the 
system.  Indeed,  as  we  will  see,  the  degree  of  Che  resultant  system  is  of 
primary  importance  in  our  theory. 


About  a  year  ago  we  gave  what  we  thought  to  be  a  correct  solution  to 
the  "frequency  domain  pole  placement"  problem  to  the  effect  that  the 


problem  admitted  a  solution  if  and  (generically)  only  if 

I 

o(q)  >  it(p)  -  1  (2.3) 


I 

3 


a 

a 


where  u(p)  is  the  total  number  of  (closed)  right  half-plane  poles  and  zeros 
of  the  plant  transfer  function,  p(s),  including  poles  and  zeros  at 
infinity.  In  the  process  of  trying  to  extend  the  result  to  the 
multivariate  case,  however,  it  became  apparent  that  2.3  was  in  error  by 
"one"  in  certain  cases.  As  such,  a  major  endeavor  during  the  past  year  has 
been  the  formulation  of  a  correct  pole  placement  theorem  for  the  single- 

variant  case. ^ 

To  formulate  the  pole  placement  theorem  we  let  p(s)  have  a  rational 
fractional  representation  as  per  equations  1.1  and  1.2  and  then  further  £ 

h 

decompose  the  numerator  and  denominator  functions  as  ffc 


n(s) 


a(s)a(s) 

m(s) 


and 


d(s) 


b(s)b( s) 
m(s) 


(2.3) 


(2. JO 


•  Here  a(s)  and  a(s)  are  polynomials  representing  the  left  and  right  half- 

plane  zeros  of  n(s),  b(s)  and  b(s)  are  polynomials  representing  the  left 
and  right  half-plane  zeros  of  d(s),  and  m(s)  is  a  Hurwitz  common 


r*- 


denominator  for  n(s)  and  d(s).  Our  pole  placement  theorem  then  results 
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from  the  following  lemma. 

LEMMA:  The  single-variate  pole  placement  problem  admits  a  solution  if 
and  only  if  the  polynomial  equation 

**  A 

a (s)x( s)  +  b(s)y(s)  =  q(s)  (2.5) 

admits  polynomial  solutions,  x(s)  and  y(s),  such  that 

o(x)  <  o(q)  +  o(a)  -  o(m)  (2.6) 

o(y)  <  o(q)  +  o(b)  -  o(m)  '.2.1) 

Moreover,  q(s)  and  r(s)  are  (polynomial)  coprime  if  and  only  if  x(s) 
and  y ( s)  are  (polynomial)  coprime. 

To  obtain  2.3  from  the  lemma  all  that  is  necessary  is  to  count 
equations  and  unknowns  in  2.5  through  2.1.  In  fact,  however,  the 
requirement  that  x(s)  and  y(s)  be  polynomial  induces  additional  constraints 
to  the  effect  that  x(s)  and  y(s)  have  non-negative  order.  Surprisingly, 
these  constraints  can  become  active  just  to  the  c-xtent  that  the  "1"  in 

equation  2.3  may  drop  out.^  Counting  equations  and  unknowns  with  this 
additional  constraint  then  yields  the  following  corrected  single-variate 
pole  placement  theorem.  liote  that  the  linear  independence  of  the  set  of 

linear  equations  defined  by  2.5  follows  from  the  coprin.enc.es  of  a(s)  and 

Ms). 

TriuJf'.LM:  The  r  ingle-variate  pole  placement  problur,  admits  a  :  elution 


and  (g-..  :u  rically)  only  if 


o(q)  ^  it(p)  -  J 


(2.8) 


where  j  is  either  1  or  0  as  indicated  in  Table  1. 


j  -  1 

j  -  o 

p(s) 

o(m)  =  o(a)  >  1 

o(m)  -  o(a)  «=  0 

and 

or 

A 

A 

proper 

o(b)  >  1 

o(b)  =  0 

o(m)  -  o(b)  >  1 

o(m)  -  o(b)  *=  0 

p(s) 

and 

or 

A 

A 

improper 

o(b)  >  1 

o(b)  =■  0 

Table  1.  Perturbation  factor  for  the  single-variate 
pole  placement  theorem. 


MULTIVARIATE  POLE  PLACEMENT:  As  a  starting  point  for  our  research  on 
the  multivariate  pole  placement  problem  we  have  considered  a  stable 
rational  multivariate  plant  characterized  by  an  n  by  n  transfer  function 
matrix,  P(s).  Since  P(s)  is  stable  it  admits  a  trivial  rational  fractional 
representation  with  numerator  P(s)  and  denominator  equal  to  the  identity. 
Of  course,  this  representation  is  trivially  (right  and  left)  coprime  via 
the  equality 

[0][P(s)]  ^  [ 1 ] [ 1 ]  =  1 


and,  as  such,  it  is  stabilized  by  compensator:  of  the  form 


(2.9) 


0(s)  =  [K(s)P(s)  +  l]_1[W(s)] 


(2.10) 


where  W(s)  is  an  arbitrary  stable  rational  matrix.  Moreover,  in  this 
special  case  the  input/output  gain  for  the  resultant  feedback  system  takes 
the  form 


H(s)  =  P(s)W(s) 


(2.11) 


In  this  highly  specialized  multivariate  case  we  formulate  the  pole 
placement  problem  by  requiring  that  W(s)  be  chosen  to  stabilize  the  system 
and  yield  an  input/output  gain  of  the  form 


H(s)  =  R(s)/q(s) 


(2.12) 


where  q(s)  is  a  prescribed  Hurwitz  polynomial  and  R(s)  is  an  arbitrary  n  by 
n  polynomial  matrix  such  that 


o(R )  <  o(q) 


(2.13) 


where  the  order  of  a  matrix  is  defined  to  be  the  maximum  of  the  orders  of 
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its  components. 

Given  P(s)  one  may  factor  out  a  (Hurwitz)  common  denominator,  b(s), 
and  then  transform  the  numerator  into  Smith  canonical  form.  After 
replacing  the  common  denominator  this  then  yields  the  fmi th-MaeMi linn  form 
for  P(s)  of  2.1  A. 


?' 


< 


**V, 


PCs)  =  U(s)D(s)V(s) 


(2.14) 


where  U(s)  and  V(s)  are  unimodular  polynomial  matrices  and  DCs)  is  a 
diagonal  matrix  with  (possibly  zero)  stable  rational  entries: 


di(s)  =  a. (s)ai(s)/bi(s) 


(2.15) 


Letting  V^(s)  denote  the  ith  row  of  V(s)  we  then  obtain  the  following 
theorem. ^ 

THEOREM:  i)  If  the  multivariate  pole  placement  problem  admits  a 
solution  then 


o(q)  >  o(bi)  -  o(ai)  -  o(Vi) 


(2.16) 


for  all  i  =  1 ,  ...  , n. 

ii)  The  multivariate  pole  placement  problem  admits  a 
solution  if 


o(q)  _>  o(bi)  -  o(a.) 


(2.17) 


for  all  i  =  1 ,  ...  , n. 

Unfortunately,  the  necessary  and  sufficient  conditions  of  2.1£  and 
2.17  do  not  coincide  except  in  the  single-variate  case  where  they  coincide 
with  the  (corrected)  single-variate  theorem. 

Needless  to  say  the  above  theorem  leaves  many  unanswered  quest ic ns. 
In  particular,  sir.ee  the  Cm ith  form  is  not  unique  the  necessary  condition 
of  the  theorem  can  conceivably  be  tightened  by  simrly  varying  the  Cr :  th 
form  employed.  Even  here,  however,  it  is  not  clear  if  a  single  Cr-.ith  i  e  r.T. 
w:  11  minimize  the  expression  of  equation  2.16  :  i:  .it ;r  1  y  ' ;  r  i 


r 

if? 

t, 


Set 


(?•- 

'j- 


> 


(per 


rally  1  doubt  it).  Furtnermcre 


» 


ip- 

I 


* 


what  happens  if  we  f; 


necessary  condition  in  terms  of  the  columns  of  U(s)? 

TIME-VARYING  SYSTEMS:  Although  many  of  the  above  techniques  can  be 
formulated  in  an  abstract  operator  theoretic  setting  their  power  and 
potential  applicability  to  the  CAD  problem  lies  with  the  simple  frequency 
domain  formulation  employed.  As  such,  during  the  past  year  we  have 
endeavored  to  formulate  a  viable  frequency  domain  theory  which  is 
applicable  to  time-varying  systems  and,  yet  still  preserves  tr  'haracter 
of  the  classical  time-invariant  frequency  domain  concept.  Of  irse,  many 
such  generalized  frequency  domain  concepts  have  been  formul;  1  over  the 
years  (including  several  by  the  principal  investigator)  using  : unctions  of 
several  variables,  matrix  and  operator  valued  functions,  etc. 

In  the  present  formulation,  motivated  by  the  work  of  Arveson  and  his 
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students,  ’  we  represent  a  time-varying  system  by  an  operator  valued 

function  of  frequency*  which,  we  believe,  has  every  analytic  property  of 

the  classical  frequency  response  for  time- invari ant  operators,  though  it 

does  not  yield  the  expected  computational  simplification  except  in  the 

.  .  .  .  .  2,15 

t .n- -invar: a he  case. 

In  its  most  general  form  (a  Hilbert  resolution  space  or  its 
generalizations)  the  required  frequency  response  takes  the  form: 


AU) 


V(u) 


-1 


«.•) 


( : .  i  s ) 


!• 

i 


r.  r.dsnt  to  r-c;  r -er-c nt  a  single  op<  rater,  the  g:v-  ti 
by  nr:  op-  r, -tor  valued  funot  ion  ( whose  v.ilue  at  0 ,  ry 
witn  t  .-'ipinal  rater-).  In  fact,  new-  ver,  wr,  t 


th‘ 


V  i 


U“d 


’  i  y  ■  r.g  :  pace  with  th.it  a  brut  the  o;  t  nr  :  n‘. 
function  which  r  :rr .i  I  tam-f-us  y  characterize:  i.ih 
tne  -.or  -u.  wr.jch  it  i:  or i in* d. 


Now,  this  yields  the  classical  frequency  response  valued  function  of  u> 
given  by 


A(tu,w)  =  [C(o)I  -  o)I  -  A)  B  +  D] 


(2.23) 


As  such,  even  though  2.23  technically  has  two  frequency  variables  ( oi  is  a 

parameter  and  oi  is  the  variable  of  a  multiplication  operator)  they  are 
clearly  redundant  and  may  be  eliminated  by  the  change  of  variable 


0)  =  0)  -  0) 


(2.2U) 


whence  we  obtain  the  classical  frequency  response  function 


A(o))  =  [C(o)I  ~  A)  B  +  D] 


(2.25) 


Although  the  above  formulation  does  not  yield  the  kind  of 
computational  simplification  associated  with  the  classical  time-invariant 
frequency  response  it  appears  to  retain  every  analytic  charact er  i  s :  i  c 
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usually  associated  with  the  classical  frequency  response.  ’  As  such,  we 
believe  that  it  can  serve  as  the  ideal  medium  with  which  to  extend  our 
control  system  design  concepts  to  the  time— varying  case  and  have  undertaken 
a  ; r- 1  in; nary  Investigation  of  this  possibility. 


A  Cl  1 V 1 : 1  h'S :  In  addition  to  the  above  described  ra  ch  we  :.nv< 


• :  it 

.  i  r.  .)•  d  .-ur  ir.v*  s  ligation  of  the  si  r  u2  tarieous  <it :  ign  problem  ’  ...nd 


initiat'd  a  ;  r-  !  is  many  invert  igatior.  of  tr.c  adept ivt  control  pret  1  or.  from 


the  point  of  view  of  cu"  design,  theory.  The  former  problem  has  proven  to 
be  extremely  hard  especially  in  the  discrete  case.  In  this  endeavor  we 
have  completed  our  abstract  algebra i s- geome  tr  i  c  char  aeter  i  za  t  i  or.  of  the 
simultaneous  stabilization  problem  which  yields  a  complete  but 
noncomputations  solution.  Furthermore,  we  have  continued  our  work  on  the 
discrete  problem,  reducing  it  to  an  (albeit  unsolvable)  interpolation 
problem.  Interestingly,  some  parallel  work  by  Byrnes,  et  al.  has  yielded 

some  generic,  but  equally  untestable  results.^ 

In  the  adaptive  control  area  we  have  formulated  an  approach  to  the 
problem  which  yields  a  solution  which  is  always  globally  stable  but  is'oaly 
locally  stable  when  the  model  for  our  plant  is  sufficiently  close  to  the 
actual  plant.  At  the  present  time  we  are  trying  to  determine  "how  close" 
is  sufficiently  close.  To  this  end  we  have  conjectured  that  the  model  must 
lie  in  the  same  connected  component  (of  an  appropriate  space  of  plants)  as 
the  actual  plant.  At  the  present  time  this  work  is  in  a  very  preliminary 
form  and,  indeed,  the  precise  definition  of  the  appropriate  "space  of 
plants"  is  still  vague. 
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